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Q\ ■ Abstract 

^ I The survival probabilities of a "prisoner" diffusing in an expanding cage 

and a "daredevil" diffusing at the edge of a receding cliff are investi- 
gated. When the diffuser reaches the boundary, he dies. For "marginal" 
boundary motion, i. e., the cage length grows as \/Ai or the cliff location 
recedes as xo(t) = —y/At and the daredevil diffuses within the domain 
x > xq, the survival probability of the diffuser exhibits non-universal 
power-law behavior, S(i) ~ t - ^, which depends on the relative rates of 
boundary and diffuser motion. Heuristic approaches are applied for the 
^ cases of "slow" and "fast" boundary motion which yield approximate ex- 

lO ! pressions for (3. An asymptotically exact analysis of these two problems 

is also performed and the approximate expressions for f3 coincide with 
the exact results for nearly entire range of possible boundary motions. 
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Consider an inebriated "prisoner" who diffuses within a one-dimensional cage defined 
by (—L(t), L(t)) and dies whenever he touches the walls (Fig. 1). We are interested in de- 
termining the probability for such a prisoner to survive until time t, S(t). Mathematically, 
the time evolution of this process in governed by the diffusion equation 

dc{x,t) _ d 2 c(x,t) 

dt dx 2 ' K ' 

within the domain —L(t) < x < L(t), subject to the initial condition c(x,t = 0) = 5(x), 
and the absorbing boundary condition c(x = ±L(t),t) = 0. Here c(x,t) is the prisoner 
density at position x and time t, and D is the diffusion coefficient. The absorbing boundary 
condition imposes the death of the prisoner whenever he touches the walls of the cage. 



1. Introduction 
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Correspondingly, the survival probability S(t) is given by the spatial integral of the prisoner 
density, 

,L(i) 

S(t) = / c(x,t)dx. (2) 

In cage of fixed size 2L, the survival probability decays as exp(— ir 2 Dt/4L 2 ) in the 
long-time limit [1]. More interesting behavior arises when we aid the prisoner by allowing 
the cage walls to recede. This obviously increases the prisoner lifetime and can also dra- 
matically change the form of the prisoner survival probability. In the following, we consider 
power-law growth of the cage length, L(t) ~ {At) a , as t — > oo. Two domains of behavior 
arise which are determined by the relative magnitudes of the basic length scales of the 
system, namely, the cage length 2L(t) and the diffusion length \l~Dt. When L(t) <C y/Dt, 
which arises when a < 1/2 or when a = 1/2 and 4<D, the prisoner diffuses faster than 
the cage walls recede and it is plausible to apply an approximation based on the assump- 
tion that the prisoner probability density is "close" to that in the case where the cage walls 
are static. This is the basis of the well-known adiabatic approximation [2]. Conversely, 
for a rapidly expanding cage, i. e., a > 1/2 ox a = 1/2 and A»D, the cage expands 
faster than the prisoner diffuses and approximations based on a "free" prisoner should be 
appropriate. 

In the interesting situation where the diffusion length of the prisoner and the cage size 
are comparable, both the adiabatic and the "free prisoner" approximations predict that 
the prisoner survival probability decays as a power law in time, S(t) ~ t - ^, but with a non- 
universal exponent which depends on A/D. To ascertain the accuracy of these heuristic 
approaches, we analyze of the diffusion equation with the moving boundary condition and 
find that the prisoner probability density may be written in terms of parabolic cylinder 
functions. From this analysis, we determine the exponent f3 for all A/D and verify that 
the results of the above two heuristic approaches are asymptotically exact. 

We also consider the related problem of an inebriated "daredevil" who diffuses in a 
one-dimensional semi-infinite domain x > xo(t), and fall to his death whenever x = xo(t) 
is reached. When the cliff location is fixed, it is well-known that the survival probability 
of the daredevil is S(t) = c(x,t) dt ~ t -1 / 2 [1]. Thus the daredevil is sure to fall off 
the cliff (although his mean lifetime is infinite). This same decay law continues to hold if 
the cliff recedes slowly, i. e., x ~ (At) a with a < 1/2 [3]. Conversely, if the cliff recedes 

2 



from the daredevil at a constant velocity v, there is finite probability for the daredevil to 
survive which rapidly approaches unity when v£q/D > 1. Here £q is the initial distance 
from the daredevil to the cliff. However, if the cliff recedes at the same rate as the daredevil 
diffuses, xo(t) ~ —\^At with A of the order of D, marginal behavior again arises in which 
the daredevil survival probability exhibits a non-universal power-law decay in time. 

It is worth mentioning that the first-passage probability of diffusion processes in the 
presence of moving, absorbing boundaries has been investigated previously by mathemati- 
cians and a substantial literature exists (see e. g., [3-5] and references therein). However, 
the methodology used in these papers is quite different from the simple-minded approach 
that we will present. In addition to differences in approach, an advantageous aspect of 
our study is that the prisoner and daredevil problems can be treated within the same 
framework. 

2. Heuristic Approaches for Prisoner Survival in an Expanding Cage 



For a fixed cage (-L, L), the solution to the diffusion equation (1) may be written as 
an eigenfunction expansion in which each eigenmode decays exponentially in time, with a 
different characteristic decay time. In the long time limit, only the most slowly decaying 
eigenmode remains and the the density approaches 

c(M) oc e-**/" 3 cos(^). (3) 

Thus the prisoner survival probability decays exponentially in time. 

Now suppose that the cage expands expands slowly, L(t) <C y/Dt. In this case, the 
adiabatic approximation shows that the density profile approaches the same form as in the 
fixed-cage case, except that the parameters in this probability distribution acquire time 
dependence to satisfy the moving boundary condition. The corresponding probability 
density is 

2L{t)) • (4) 
with the amplitude f(t) to be determined. Substituting Eq. (4) into Eq. (1) leads to 

• f Dn 2 \ „ / 7TX \ . TlX . • 



When L(t) grows as (At) a with a < 1/2, the second term on the right-hand side may be 
neglected and the leading behavior of the amplitude f(t) is given by 



/(£) -> exp 



n 2 D 



dt'L- 2 (t') 
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(6) 



The full asymptotic behavior of f(t) is expected to also contain a power-law prefactor in L. 
However, this prefactor is not accessible within our naive approach. The leading behavior 
of fit) now gives 

S{t)= [ m c(x t t)dx~ - f(t)L(t). (7) 

J-L(t) * 

Thus for a < 1/2 the leading behavior of the survival probability decays as a stretched 
exponential in time 



S(t) -> exp 



(8) 



4(1 - 2a)A 2a 

For the marginal case of a = 1/2, the second term in Eq. (5) is no longer negligible. 
Following the above procedure nevertheless, we find a non-universal power-law behavior, 

S(t)~t-e, with 0=^. (9) 

Given the nature of the approximation employed, this prediction is anticipated to be valid 
in the limit i<D. 

In the complementary case of a rapidly growing cage, L(t) 3> y/~Dt, i. e., a > 1/2 or 
a = 1/2 and A 3> D, it is plausible to assume that the prisoner density profile approaches 
a Gaussian but with a decaying overall integral. Thus we are led to hypothesize 

c (x,t) ~ ^SLexp (-^-] . (10) 

Although this distribution does not satisfy the absorbing boundary condition, the inconsis- 
tency is expected to be negligible, since the density is exponentially small at the cage walls. 
The decay of the mass may be found by equating the flux to the cage walls, 2j = —2D^, 
with the mass loss. An elementary computation shows that the survival probability ap- 
proaches a constant for a > 1/2. In the marginal case of a = 1/2, 



S = -?iAexp 4). (11) 



t V 4nD V 4L> 



This again leads to power law behavior for the survival probability, S ~ t 13 with 

" = < 12 > 

The predictions of this "free prisoner" approximation should be accurate when A»D. As 
is intuitively expected, when A/D decreases toward 0, /? diverges. This corresponds to the 
survival probability exhibiting faster than power-law decay, as predicted by the adiabatic 
approximation. Conversely, as A/D — > oo the prisoner diffuses more slowing than the cage 
expands and there is a finite chance to survive asymptotically, i. e., f3 — > 0. 

3. Asymptotic Analysis for the Marginally Growing Cage 



Let us now investigate more carefully the borderline case where the cage grows as 
Lit) ~ \J~A±. Within a scaling formulation, it is natural to hypothesize that the density 
can be written in terms of the dimensionless variables 



x x 



as t — > oo. Since both the cage length and the diffusion length grow at the same rate, it 
proves convenient to consider the basic variables to be £ and p = £/cr = A/D and write 
the concentration as 

c(M) ~ f _/, - 1/2 Cp(0- (13) 

The power law prefactor is chosen to ensure that the survival probability decays as t~@, 
as defined previously. 

Substituting Eq. (13) into Eq. (1), we find that the scaling function C p (£) satisfies the 
ordinary differential equation 

Dd 2 C \AC ( a \\ n n 



a <ie 2 s d£ 

where the p dependence has been dropped, for notational simplicity. Introducing £ = 
r] \J1D/A and C(£) = e~ v / 4 V(rj), transforms Eq. (14) into the canonical form for the 
parabolic cylinder equation [6] 



A spatially symmetric solution to Eq. (15) (appropriate for the prisoner starting in the 
middle of the cage) is 

V{r 1 ) = \{V w {r 1 ) + V w {-r 1 )), (16) 

with V u {ji) the parabolic cylinder function of order v. Finally, the relation between the 
decay exponent (3 and A/D is determined by the absorbing boundary condition 




= 0. (17) 

One can easily determine the limiting behaviors of (3 = f3(A/D) for A/D — > and 
A/D — > oo and thus check the validity of the heuristic predictions given in Eqs. (9) and 
(12), respectively. In the former case, the exponent (3 is large and hence, from Eq. (15), 
the profile approaches the cosine form. Thus Eq. (9) provides the correct asymptotics. In 
the latter case, A ^> D, (3 — > and Eq. (15) approaches the Schrodinger equation for the 
ground state of the harmonic oscillator, for which T>o(r/) = exp(— 1] 2 /4). For small but 
finite (3 it is natural to seek a perturbative solution 

V(rj) = exp(-77 2 /4) + /M(??) + . . . . (18) 

Substituting this expansion into Eq. (15) yields an inhomogeneous linear equation for the 
correction A(rj): 

Introducing B(rj) through A(rj) = exp(— rf / A)B{rf) , we obtain B" — rjB' = —2. By solving 
this latter equation, the perturbative solution for T>(rj) is 



Jo Jo 



(20) 



Combining Eq. (20) and the absorbing boundary condition V(^J A/2D) = we reproduce 
(after a straightforward but lengthy computation) the asymptotics given by Eq. (12). 
Thus we have rigorously justified the previous heuristic predictions predictions for the 
decay exponent (3. 

It is also interesting to consider the mean lifetime of the prisoner, (t) = J °° dtS(t). 
This quantity is finite for (3 > 1 and infinite for (3 < 1. The borderline case of (3 = 1, 
corresponds to the second excited state of the wavefunction in the harmonic oscillator 
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potential in Eq. (15), for which the solution for V(rj) is V{rf) = (1 — i] 2 )e~ ri / 4 . The 
boundary condition of Eq. (17) now gives A = 2D. Thus the borderline case between a 
finite and an infinite mean survival time corresponds to A = 2D. It is gratifying that for 
this case of A = 2D, the simple-minded adiabatic approach gives f3 = 1.234. This provides 
a sense for the accuracy of the adiabatic approach in the regime A < 2D. 

4. Survival of a Daredevil at the Edge of a Cliff 



Let us now turn to the case of an inebriated daredevil who survives if he remains 
within the semi-infinite domain xo(t) < x < oo. As mentioned in the introduction, the 
most interesting case is where xo(t) ~ — \J~M as t — > oo, corresponding to the cliff receding 
at the same rate at which diffusion tends to transport the daredevil to the cliff. This 
situation can be analyzed by methods similar to those applied for the prisoner in the 
marginally growing cage. 

It is convenient to first change variables from (x,t) to (x' = x — xo(t),t) to fix the 
absorbing boundary at the origin. Thus the initial diffusion equation is transformed to the 
convection-diffusion equation (where the prime is now dropped) 

dc in dc ^ d 2 c „ , , 

with the fixed absorbing boundary condition c(x = 0, t) = 0. As in the case of the prisoner 
in the marginally expanding cage, we apply the same scaling assumption, Eq. (13), for the 
probability density of the daredevil. Substituting this form into Eq. (21) gives 

Dd 2 C 1_ ^dC ( n 1\„ 

+ o(£-l)37+ /*+- C = 0. (22) 



A d£ 2 2 ' d£ V 2 

Transforming from ^ = — ^ and C(^) to r\ and T>(rj) in a slightly different way than 
previously, 

e- 1 = V^^ C(0 = exp(-^P(r 7 ), (23) 

we find that T>(rj) satisfies the same parabolic cylinder equation (Eq. (16)) as in the prisoner 
problem. However, slightly different boundary conditions apply. Death of the daredevil at 
the edge of the cliff implies 

V(-y/A/2D) = 0. (24a) 



On the other hand, S(t) = J °° dx c(x,t) < 1 implies the boundary condition at 77 = oo, 

V{r} = 00) = 0. (246) 

Mathematically, the determination of (3 and T>(rj) is equivalent to finding the ground 
state energy and wavefunction of a quantum particle in a potential composed of an infinite 
barrier at 77 = — ^/A/2D and the harmonic oscillator potential for 77 > — y / A/2D. Higher 
excited states do not contribute in the long time limit. A general solution to Eq. (22) 
satisfying V(oo) = is 

V( V ) = V 2f3 (rj), (25) 

and the absorbing boundary condition T>2p{— ^ A/2D) = determines the relation be- 
tween the decay exponent (3 and A/D. 

Since Eq. (25) provides only an implicit relation (3 = f3(A/D), it is useful to determine 
the limiting behaviors of (3 for small and large values of A/D. These relations can be 
derived directly from Eq. (21) and elementary facts about the quantum mechanics of 
the harmonic oscillator, rather than relying on mathematical properties of the parabolic 
cylinder functions. In the limit of slowly moving boundary, A < D, the wall is close to 
the origin in the 77-variable. When the wall is exactly at the origin, the ground state of the 
truncated potential is obviously the first excited state for pure harmonic oscillator, namely, 
2/3 = 1 and V(rj) = 77 exp(-?7 2 /A). For A < D, we therefore expect that < 1 - 2j3 < 1. 
This again suggests the perturbative solution, 

V(rj) = 77 exp(-77 2 /4) + (1 - 2j3)A{rf) + .... (26) 

Substituting this expansion into the parabolic cylinder equation yields for .4.(77): 

v + G-ty-** H - (27) 

Introducing B(t]) through A(rj) = 77 exp(— ?7 2 /4) £>(r?) we find B" — (77 — 2/r])B' = 1. Solving 
this equation subject to the boundary condition Eq. (24b) gives B(rj), from which one 
ultimately obtains 



V{rj) = 77 e" 7 ' 2 / 4 + (1 - 2/3) 77 e"" 2 / 4 / d m 7]^ 2 e^ /2 j dr] 2 rfi ^ l ' 2 + ■■■■ 



(28) 



Applying the absorbing boundary condition Eq. (24a) gives 

— . (29) 
1 2 V AtxD v ; 

One can treat the opposite limit A ^> D similarly. In terms of the coordinate rj, the 
location of the wall goes to — oo. Hence the unperturbed ground state for this system is just 
the ground state for the pure harmonic oscillator, namely, (3 = and V{rf) = exp(— if /4). 
Following the same perturbative approach as in the complementary case of i C D, we 
find 

/•OO /-OO 

V{rj) =e- r > 2/4 + 2/3e- r > 2/4 / d^e^ 2 dr) 2 e~^ /2 + 0(/3 2 ). (30) 

Combining Eq. (30) with the absorbing boundary condition one finds the same expression 
for (3, Eq. (12), as was found for a finite cage. 

5. Summary and Discussion 



We have presented a heuristic and an asymptotically exact approach to determine 
the survival probability and the density distribution for (i) a prisoner in a growing cage, 
(— L(t), L(t)), and (ii) a daredevil in the domain x > xo(t) with a cliff at x = xo{t). We 
were primarily concerned with the "marginal" case where L(t) = \f~At and x$ = —\f~At 
(with A of the order of D), so that boundary of the system recedes at the same rate at 
which diffusion tends to bring the diffuser (prisoner or daredevil) toward his demise. In 
these marginal situations, the survival probability of the diffuser exhibits a non-universal 
power-law decay in time. The value of the decay exponent in the limiting cases of A 3> D 
and i € D can be obtained by simple arguments. These limiting behaviors are found to 
coincide with the results from an asymptotic analysis of the underlying equation of motion. 

For the prisoner problem, our results can be straightforwardly extended to general 
spatial dimension. Following the same adiabatic and "free prisoner" approximations that 
were applied in one dimension, we find that the decay exponent (3 becomes 

n 2 D 

„ l&T- adiabatic; 

P= 1 , A, ^ ( 31 ) 



r(f) 



w) e M~w) free - 



Here jd is the first positive root of the spherical Bessel function Jd/2-i(x)- Similarly, in 
the case of marginal cage growth, the (/-dimensional analog of the scaling ansatz, Eq. (13), 



can be applied, leading to a generalization of Eq. (14). In three dimensions, in particular, 
the additional transformation C(£) = leads to the same parabolic cylinder equation 

(14) for JF, but with the parameter (3 + \ replaced by (3 + 1. 

It is instructive to compare the behavior of the survival probabilities given here with 
those of the related problem where the absorbing boundaries themselves undergo diffusive 
motion. These are situations for which exact solutions have been given previously [7]. 
For example, consider the survival of a diffusing daredevil when the position of the cliff 
also diffuses with a diffusivity A. This situation is trivially isomorphic to the case of a 
static cliff and a daredevil with diffusivity D + A. Thus the survival probability decays 
universally as t -1 / 2 . As might be expected, a cliff which is systematically receding from a 
diffusing daredevil with L(t) oc \ft leads to a larger survival probability compared to the 
case of a stochastically moving cliff. 

On the other hand, the survival of a diffusing prisoner inside a cage where both walls 
diffuse (each with diffusivity A) is more interesting. A variety of exact solutions show that 
the prisoner survival probability decays non-universally, S(t) ~ t _/3 , with decay exponent 
(3 = n/2cos~ 1 (D/(D + A)) [7]. For rapidly diffusing walls, (3 — > 1, while the corresponding 
limit for systematically receding walls is (3 — > 0. Clearly if cage walls are receding rapidly, 
the prisoner is more likely to survive compared to the case where the cage walls are diffusing 
rapidly. On the other hand, in the limit of slowly diffusing walls, (3 — > \J Strangely, 
this is almost the square-root of the corresponding expression for (3 quoted in Eq. (9). 
Intriguingly, the prisoner is more likely to survive in a stochastically and slowly growing 
cage than in a cage which grows systematically at the same average rate. 
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Figure Caption 

Figure 1. (a) The prisoner in the expanding cage, and (b) the daredevil at the edge of a receding 
cliff. 
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